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A zero-field Ising model with ferromagnetic coupling constants on the so-called Labyrinth 
tiling is investigated. Alternatively, this can be regarded as an Ising model on a square 
lattice with a quasi-periodic distribution of up to eight different coupling constants. The 
^ I duality transformation on this tiling is considered and the self-dual couplings are determined. 

^ I Furthermore, we analyze the subclass of exactly solvable models in detail parametrizing the 

O ! coupling constants in terms of four rapidity parameters. For those, the self-dual couplings 

^ I correspond to the critical points which, as expected, belong to the Onsager universality 

Q\ ' class 

1 Introduction 

O ■ The understanding of phase transitions in 2D systems has considerably increased during 
c/2 ■ 

the last decade. In particular, a lot is known about the conformal structure at the critical 



point of a whole hierarchy of models [|T|. On the other hand, numerous series of solvable 
vertex models 0, |^ and IRF (Interaction- Round-a- Face) models 0, ^ have been constructed. 
Furthermore, the so-called chiral Potts models [^ ^ attracted a lot of interest, mainly due to 
the existence of superintegrable cases |^ (see also [0 and references therein) and applications 
to 3D systems p|. 

While many of these models can be solved in the presence of a temperature-like deviation 
from criticality, solvability is in general restricted to the zero-field case (apart from e.g. the 
hard hexagon model p| which is a limiting case of the triangular Ising model in a symmetry- 
breaking field) . Recently, a series of models which can be solved in the presence of a symmetry 
breaking "field" was found [Q , among them a model which belongs to the universality class 
of the Ising model in a magnetic field. But even on regular lattices, notably the square 
lattice, the word "solvable" does not mean that all physically interesting quantities have been 
calculated explicitly, in particular the computation of two-point (or even higher) correlation 
functions still appears to be a hopeless task in the general case. 
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Even less, however, is known about models on non-periodic graphs. The main reason 
is that quite a number of new phenomena show up: already the Lee- Yang zeros of the 
classical ID Ising model on substitution chains reveal gap structures that are unknown in 
the periodic case |10|. Similarly, the Ising quantum chain on such substitution structures 
n and the 2D Ising model with aperiodicity in one direction W^, K^ resemble 



[0,0 

the periodic situation at its critical point only under special conditions. These examples 

also show that results from conformal field theory cannot generally be applied to systems 

on non-regular graphs which is not too surprising as they rely on the basic assumption of 

conformal invariance of the system at the critical point. 

The restricted knowledge about non-periodic systems is due to the significantly increased 

difficulty to reach rigorous results, the study of 8-vertex models on relatively general grids 

[p^ , |18[ or certain IRF models on (rhombic) duals of regular de Bruijn grids |jl8|, |T9| being 

untypical exceptions (see |2^ for a definition of the grid formalism). Other investigations 

make use of various approximative techniques such as, for instance, real space renormaliza- 

tion group ||2I|, ^, mean field calculations [^, computer simulations |^ |2S|, |2^, PB| , 

series expansions ||29|, and finite-size scaling methods ||3^, but almost exclusively deal with 



Ising or Potts models on Penrose tilings. 

In this article, we consider the zero-field Ising model defined on a different quasiperiodic 
graph in 2D. We have chosen the so-called Labyrinth tiling, see |^ for a detailed geometric 
description. It turns out that this tiling is self-dual as a graph. This allows - under certain 
assumptions such as uniqueness of the critical point - a simple argument to obtain infor- 
mation about the location of critical points, though it turns out to be incomplete. In the 
case of the frequently studied Penrose patterns, duality and correlation inequalities are less 
conclusive, but have been used to obtain bounds for the critical temperature P^ , 

From the Peierls argument 



5^, ^ it is clear that there must exist at least one phase 
transition in our system. We do not expect that the continuous transition in the Ising model 
on regular lattices becomes a first-order transition in our case since introducing "randomness" 
generally tends to smooth out transitions, though it is certainly not completely obvious (e.g., 
the four-spin interaction Ising model on the tetrahedra of the 3D fee lattice in fact has a 
first-order transition at its self-dual point [^ ). 



Unfortunately, we could not decide the question of uniqueness of the critical point in 
the ferromagnetic case in general, though one might expect this to be true if all coupling 
constants are strictly positive and uniformly bounded from below. To further attack this 
problem, we have investigated the solvable cases of the model by means of the rapidity line 
approach [^. We find a 4D manifold in the 8D space of couplings where the Ising model 
on the Labyrinth is solvable. Therein, one finds a 3D critical surface, hence it contains one 
temperature-like parameter that describes the deviation from criticality. The result is given 
explicitly by means of some geometric properties of the Labyrinth. In all these cases, the 
critical point is unique and belongs to the Onsager universality class as expected [0, ^ . 



2 The Labyrinth tiling 

To investigate exactly some typical properties of generic non-periodic graphs in 2D, we need 
a simple example which is more complicated than the Cartesian product of a Fibonacci chain 
and a ID periodic lattice (which is what underlies the well-studied Ising quantum chains), 
but still considerably simpler than a Penrose tiling. The natural choice is the so-called 
Labyrinth tiling ^1|], which got its name from the properties of certain colourings. There 



exist several ways to construct it, see 0, 39 1 for a detailed description. For our purpose 



an approach via the Cartesian product of a ID substitution sequence with itself is most 
appropriate which we will now briefly outline. 

2.1 The silver mean substitution rule 

Let us start with the generation of a ID chain, the so-called silver mean chain, by repeated 
apphcation of the two-letter substitution rule 

''■ h ^ hah ^'•') 

to the letter a. The corresponding substitution matrix R is obtained from the standard 
abelianization process, compare [^ and refs. therein, and reads 




R=\ . o • (2-2) 



The Perron- Frobenius eigenvalue is the silver mean, A = 1 + -\/2, the corresponding eigenvec- 
tor reads (1, A)*, or ((3 — A)/2, (A — l)/2)* in statistical normalization. Since R is symmetric, 
left and right eigenvectors are the transpose of each other. Thus, the entries have two mean- 
ings: on the one hand, they tell us the relative frequencies of a's and 6's in the infinite chain 
- the ratio of which is the irrational number A wherefore the chain cannot be periodic. On 
the other hand, they give us the natural choice for the lengths of two intervals with which 
we can represent the chain geometrically such that the substitution rule ( |2.1| ) gives rise to a 
inflation/deflation symmetry. In this setup (which we shall use for graphical presentation), 
a stands for a short interval, while that attached to 6 is a factor of A longer. 

2.2 Construction of the Labyrinth 

The next step is to generate a 2D tiling. First, we take an orthogonal Cartesian product 
of two identical silver mean chains in the proper geometric representation. This way, we 
obtain a quadrant with an orthogonal grid, but two different spacings in each direction. 
Now, starting from the lower left corner, we mark every second vertex point of this grid and 
connect them with the nearest neighbours of the same kind. This would then give a fourfold 



symmetric, non-periodic tiling of the plane (after filling the 3 free quadrants with rotated 
copies of the grid constructed) . 

In what follows, we shall, however, mainly need so-called periodic approximants. They 
can easily be obtained from the periodic approximants of the silver mean chain. The n-th 
iteration of the substitution rule (remember that we start with the letter a for n = 0) gives 
a finite chain with g^ cells with the generalized Fibonacci numbers g^ = 1, gi = 1, and 



gn+i = '^Qn + Qn-ii comparc eqs. (4.19) and (4.20) in Ref. |^. Now, due to the special 



choice of the substitution (|2.1| ), we can close the chain periodically not only after the full 
number of cells, but also after gn — ^ cells, which is always an even number. 

If we do that for our finite grid (obtained from the Cartesian product of two identical finite 
chains), we wrap it on the torus and get the periodic approximants we need, compare Fig. 1. 
They have the nice property that neither new vertex configurations nor other mismatches 
occur. Such an approximant, obtained from the n-th iteration step of the substitution rule, 
is called £„. Note that the number iV„ of vertices in £„ is given by 

Nn = lign-iy. (2.3) 

The tiling was constructed from half the points of the orthogonal grid. If we use the 
other half, we obtain another tiling of the same kind, and, even more, the dual graph to the 
Labyrinth tiling. This also establishes the checkerboard structure we will need later. The 
duality of the two graphs remains true for the periodic approximants, and Fig. 2 shows both 
C^ and £* for n = 4. It is also obvious that the Labyrinth tiling is topologically a square 
tiling - but we have now performed a non-periodic deformation. Furthermore, this tiling 
is quasiperiodic 131], p9l, but it can be seen as a modulated structure and is therefore not 



a true example of a quasicrystal [41]. But this does not matter for our purposes as many 
arguments, in particular the rapidity line approach of Sec. 4, apply also to discrete, locally 
finite structures obtained from regular grids via dualization, compare |1^, [1^, 0. Note, 



however, that these examples are usually not self-dual as a graph. 

Before we proceed, let us describe some properties of the Labyrinth tiling which are useful 
in our present context, for additional material we refer to [^, |3^. First of all, the tiling 
is built from three tiles, a square (abbrev. by the letter A from now on), a kite (B) and a 
trapezoid (C). They show altogether three different edge lengths which we call i, m, and s 
for long, medium, and short, respectively. We will now summarize their statistics. 

2.3 Tiles and vertices 



The substitution rule (|2.1| ) of the underlying ID chain induces a substitution rule for the 2D 
tiles which is local, compare Fig. 3. The combinatorial net rule reads 



:2.4) 



A - 


-^ A + AB + AC 


B - 


-^ A + 3B + 2C 


C - 


-^ A + 2B + C. 



This rule can be summarized by the tile substitution matrix Rt 



Rt 



(I 1 1\ 

4 3 2 

4 2 1 



:2.5) 



compare fig. 2(b) of [^. Since this matrix is no longer symmetric, one has to distinguish 
geometric and statistical quantities. We have chosen the convention where the statistical 
properties can be read from the usual (right) eigenvectors. 

The Perron-Frobenius eigenvalue is A^ = 3+2a/2 with eigenvector (5— 2A, 6A— 14, 10— 4A)* 
which shows the frequencies of the tiles in the thermodynamic limit. Due to self-duality of 
the tiling as graph, this tells us also the frequencies of the vertex configurations, see Fig. 
4. Here, the four orientations of the kite and the trapezoid (resp. that of the corresponding 
vertex configurations) are summed over, respectively. Since the tiling has fourfold symmetry 



31|], each orientation is equally frequent which we will use later. 



2.4 Statistics of bonds 

For the Ising model to be studied next, we need some information about the bonds of the 
tiling. One important observation, compare Figs. 1 and 2, is that transversally intersecting 
dual bonds of C and C* are always of the same type (i.e., d intersects with t etc.). Further- 
more, this duality is still consistent with an anisotropic generalization, where we distinguish 
altogether eight different types of bonds. This holds true both for the thermodynamic limit 
and for the finite approximants. 

Also, we can calculate the frequency of the bonds. The substitution rule for the tiles, 
summarized in ([^.5|), induces the following edge replacement rule. 



From this, one can read m 
tution matrix i?^ reads 



a. + m - 

i + 2m + s - 
2i + m and s 



U + Am + s 
6i + 5m + s 
9i + Qm + s. 



(2.6) 



wherefore the corresponding bond substi- 



Rh 



(A 

4 



1\ 





(2.7) 



Here, the Perron-Frobenius eigenvalue is again A, now with eigenvector (|, A — 2, | — A)* in 
statistical normalization, i.e., exactly 50% of the bonds in the thermodynamic limit are of 
type i. etc. Again, different orientations of the same type of bond are summed over, but are 
equally frequent due to the symmetry of the tiling. We have now all geometric data we need 
for the description of the Ising model. 



3 Ising model and duality 

We now consider a zero-field Ising model with spins a G {1,-1} located on the vertices 
of the periodic approximant £„. Let J^, J^, Js denote ferromagnetic coupling constants 
for neighbouring spins connected by a long, medium, or short bond, respectively. In other 
words, a bond of type x between neighbouring spins at locations i and j contributes —Jx<Ji<Jj 
to the total energy, where Jx > and x G {i, m, s}. 

The canonical partition function Z„ of the periodic approximant £„ is the sum over all 
configurations a on C^ and has the form 

^n = Y.\ H exp{K^aiaj) [] ^^M^m^i^^j) IT exp{K^aiaj) . (3.1) 

o- \<i,j>i <ij>m <i,j>s J 

Here, <i,j >x denotes all nearest neighbour pairs at positions i and j which are connected by 
a bond of type x and K^ = pjx, x G {i, m, s}, where P = l/ksT is the inverse temperature. 
In complete analogy to the treatment of the square-lattice Ising model in [^ (we use the 
same notation throughout the argument) we can write down both low- and high-temperature 
expansions for the partition function ( p.lD in terms of polygons on the lattice £„. Let us 
commence with the high-temperature series. 

3.1 High-temperature expansion 

Consider an arbitrary, but fixed £„. Since ctjCTj G {—1,1}, one has 

exp(i^o"jCrj) = cosh(_ft') (1 -|- tanh(i^)crjcrj) . (3.2) 

Now, we introduce, for x G {i, m, s}, 

Vx = tanh(ir^) (3.3) 



insert this and the previous expression into ( |3.1| ), and expand the products. If rux denotes 
the number of bonds of type x in C^, we obtain 2™*+™'"+™^ terms of the form 



N„ 
-"m '^ s 



v.^vlTv? ■ n <^ , (3.4) 



i=l 



where Tx is the total number of x-bonds in the polygonal representation of the corresponding 
configuration, compare |@|, and z/j is the number of lines with site i as an endpoint. Since 
ai G {-|-1, —1}, summation of (|3l^) over the spin configurations will not contribute to Z„ 



unless all Vi are even, when it sums up to vYv^vl" ■ 2^". 

This allows us to rewrite the partition function (|3.1|) as follows 

Zr, = 2^" (cosh ir^)"^ (cosh i^^)"'" (cosh i^J*"^ J2 ^7^m<' (3-5) 

PcC„ 



where the summation is now performed over all polygon configurations P on £„. Note that 
(|3.5|) is the exact expression for any (finite) periodic approximant because we do not truncate 
the series - the name only stems from the importance ordering w.r.t. high temperature. We 
will now turn to the appropriate counterpart for low temperature. 

3.2 Low-temperature expansion 

Let nix be as above and denote by r^ the number of unlike nearest-neighbour spins linked 
by a bond of type x. Then, for a given configuration, there remain me — ri spin pairs of type 
i etc., wherefore the corresponding term in Z„ takes the value 

exp(ir^(m£ - 2re) + K^{mm - 2r^) + K^ms - 2r,)) . (3.6) 

Now, if any two adjacent spins in C^ are different (and only then), we draw the dual bond 
in £* which is of the same type. This gives us a set of precisely ri,rm,rs lines (bonds) of 
type £,m,s in £*, respectively, which form closed polygons. 

Since, in turn, each such polygon configuration represents exactly two spin configurations, 
the partition sum (|3.1|) can be rewritten as M 



Zn = 2 exp{Kf,m(, + K^mm + K^nis) J2 exp (-2/C^r^ - 2ii'„r„ - 2ir,rs) . (3.7) 

The sum is now over all polygon configurations P in the dual lattice £*, which, however, is 
isomorphic with the lattice £„ itself wherefore the sum over £* can be replaced by one over 
£„. Then, we can directly compare with Eq. ( |3.5| ) above because the meaning oim£,ri etc. 
is now the same. 

Note that there is a slight difference from the square-lattice model in where the 
horizontal and vertical lines interchange by going to the dual lattice. In our case above, we 
only associated different coupling constants to bonds of different length. We will discuss an 
anisotropic generalization with eight different (ferromagnetic) coupling constants later. 

3.3 Free energy 

The average free energy per site /, respectively the corresponding dimensionless quantity ip, 

i^f = ^ = -il- ]^ln^n (3.8) 
can now be expressed in two ways using the above expansions ( |3.5| ) and ( |3.7|) 

-V = K,fie + K^fim + K^fis + $(e-2^^e-2^™,e-2^0 (3-9) 

= ln[2(coshir,)^Mcoshirj'^-(coshirj'^^] + ^v,,v^,v,) . (3.10) 



Here, we introduced the notation 

/i. = „lim^p (a;G{£,m,.}) (3.11) 

and 

1 ( \ 

*K^„^,^s) = i™o^^^ ^ <'^m™<M (3-12) 

Thus /Xa; is the frequency of the bonds of type x relative to the number of vertices, which 
means that the /i^,. are normahzed according to 

fil + fim + fJ-s = 2 . (3.13) 

This follows from the co-ordination number of any vertex being Uc = 4 and each bond 
belonging to precisely two vertices. The actual values follow from the remarks after Eq. ( p.7|) 
and are given by 

fie = I, fim = 2X-4, /x, = 5 - 2A (3.14) 

with A = 1 + 72. 

3.4 Duality and critical point 

We introduce dual couplings K* by 

exp(-2irj = tanh(ir:) = < (3.15) 

(or, equivalently, by exp(— 2K*) = tanh(ii'^)) with x G {i,m, s}. They satisfy 

sinh(2ii'^) sinh(2ir*) = 1 (3.16) 

which is a more symmetric way to present the relationship between the two sets of couplings. 
Also, Eq. ( |3.16| ) shows directly that the duality transformation 

* : K, ^ K: (3.17) 



is an involution on IR^, where it is also analytic. Using Eq. ( |3.13D , one obtains from Eqs. 
(|3.9|) and ( |3.10|) the following transformation of ip under *, 



^(i^;,ir;,i^:) = ^(ir„ir„,irj + lln[(sinh2ir,)'^^(sinh2irj'^-(sinh2irj'^^] . (3.18) 

Therefore, the duality transformation relates the free energies at different couplings which 
implicitly depend on temperature. In general, they do not belong to the same Ising system, 
wherefore we cannot say much about phase transition points (= non-analyticity points of 
ip). However, if we define 

M := {{Ke, K^, K,) e Rl \ ^{Ke, K^, K,) not analytic } (3.19) 



we know that Ai is invariant under duality transformation *. It would now be interesting 
to know the hierarchy of ^-invariant subsets of A4 because this would tell us the critical 
structure of the model, but, unfortunately, we do not have enough information about A4 to 
do so. Let us therefore at least look for invariant points. The only fixed points of Eq. ( |3.17 ) 
are obtained from 

sinh(2is:^) = sinh(2is:„) = sinh(2irj = 1 , (3.20) 

i.e., Kf^ = K^ = Kg = ln(H-A/2)/2 = 0.44068679 . . .. This, of course, is just the phase tran- 
sition point of the periodic square lattice Ising model ^, If. Unfortunately, this does not 
tell us anything about the critical behaviour in the non-periodic case. 

Before we continue, let us remark that the additional term on the right-hand side of 
Eq. ( p.l8| ), which is anti-symmetric under duality transformation, vanishes at the self-dual 



point. Therefore, the latter lies on the surface defined by 

(sinh2i\:^)^*(sinh2i^^)^'"(sinh2irj^^ = 1 (3.21) 

with /i^, /im, and //« as listed in Eq. ( p.l4| ). This surface consists of all points in coupling space 
for which the free energy is invariant under duality transformation (though the couplings 
may change). It would be interesting to see how phase transition points for non-periodic 
models (which certainly exist for non-vanishing couplings due to the Peierls argument, see 
3^], ch. V of [Q and, for the proper generalization needed here, ch. V.5 of IQ) are located 



relative to this surface, i.e., how M. is related to it. 

As mentioned before, we do not know the full solution for the non-periodic model. In the 
ferromagnetic regime, one might expect a single critical point of the Onsager universality 
class, but even this appears to be quite complicated to decide. 

If one also allows anti-ferromagnetic interactions, the general picture will, of course, be 
different due to the presence of frustration, compare [Q (see also |3^, for the effects of 
frustration even in the ID case). Things get also more complicated then by local energetic 
degeneracies, a whole hierarchy of which can be seen as a consequence of inflation/deflation 



symmetries [^6|, ^. However, if a// interactions are anti- ferromagnetic, the behaviour of the 
system will be the same as in the purely ferromagnetic case since the Labyrinth is a bipartite 
graph and an overall sign can be absorbed reversing the spins on one sublattice. 

3.5 Anisotropic generalization 

So far, we have dealt with three types of bonds only. Like in the square lattice case, compare 
[0 , we can however distinguish "raising" bonds from "lowering" bonds. If we introduce the 
label xy for the box with abscissa x and ordinate y and K^y = pj^y {L^y = pj~y) for the 
couplings of a raising (lowering) bond in it (see Fig. 5), it is easy to derive from Fig. 2 that 
the proper duality transformation^ is now given by 

^ ■ I ^y I I ^ I ^y \ (3.22) 

\ ^xy J \ Lly j 



^We will use the same symbol since misunderstandings are unlikely. 



with K*y and L^^ defined by 

sinh(2K*^) sinh(2L,^) = 1 

sinh(2L:^) sinh(2ir,^) = 1 . ^ ' ' 

Here, xy is any element of tfie set of labels {aa, ab, ba, bb} where a and b refer to ( p.l| ) and 
hence to a short and a long interval, respectively, of the original ID chain, see Fig. 1. Again, 
* is an analytic involution, this time on IR^. 

Also, one can use the definition ( ^.191 ) w.r.t. this 8D space. From the Peierls argument. 



one would expect the set M. to be some sort of codimension one manifold, possibly with a 
rather complicated structure, where again the knowledge of ^-invariant subsets is desirable. 
Let us look at the submanifold S of self-dual points, which is easily calculable. With the 
abbreviation 

S^y := sinh(2ir,^) sinh(2L,^) (3.24) 

one obtains the fixed points of Eq. ( p.22[ ) as the solutions of 

Saa = Sab = S^a = S^b = 1 • (3.25) 

A little later, we will relate this self-duality surface S to critical points obtained from exactly 
solvable cases. 

In order to identify the 4D manifold S defined by Eq. ( |3.25D with a surface of critical 
points, one needs the following arguments and assumptions. For every point ^ G iS one can 
easily find a curve Kxyij), L^yij) parametrized by r G IR+ with the following properties: 
(i) Kxyir) and L^yij) are strictly increasing functions of r, (ii) the curve is mapped onto 
itself under duality, (iii) for some (unique) tq, the curve passes through ^, (iv) K^yiS^) = 
Lxy{0) = and Kxy{oo) = L^y^oo) = oo. Obviously, such a curve need not belong to a 
single Ising modelQ. Nevertheless, r plays the role of an inverse temperature wherefore the 
usual arguments apply: (a) there is at least one phase transition along the curve, compare 
[P^] , and (b) assuming uniqueness, the phase transition point coincides with the self-dual 
point ^. We shall say more about criticality in Sec. 4. 

Of course, one can also work, step by step, through the generalizations of the equations 
of the previous sections to find the transformation of the free energy under * . The manifold 
of couplings where the free energy is invariant under duality is now given by 

iSbbr^HSbar-^HSabr-^HSaar-^' = 1 (3.26) 



which contains iS as a submanifold. In Eq. ( |3.26| ), we have used the frequencies of the bonds 
and their isotropic distribution over the different orientations. 

Again, it would be instructive to see how the various surfaces are interrelated, which 
we will now partially explore for a subspace of coupling space where the model is exactly 
solvable. Here, the set of critical points in fact coincides with the corresponding submanifold 
of the self-duality surface. 

^Remember that this is defined through a fixed set of coupling constants J^ . 
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4 Rapidity lines and solvability 

The checkerboard structure of the original graph gives us the opportunity to introduce 
so-called rapidity lines for the parametrization of the couplings, compare [^ for details. 
Whereas Baxter considered the most general case of a "Z-invariant" Ising model [^ (which 
is a special case of the exactly solvable "Z-invariant" zero- field eight- vertex model Ijl^), we 
are only interested in models which have at most eight different coupling constants associated 
with the eight different bonds in the Labyrinth tiling. This leads to exactly solvable models 
with commuting transfer matrices and couplings K^y and Lxy parametrized by differences of 
four rapidity parameters u^ and u' (where both x and y can be either a or h) for the horizontal 
and vertical rapidity lines, respectively, and one temperature-like variable Vt which describes 
the deviation from criticality. Altogether, this means that we have a 4D solvable subspace of 
the 8D space of couplings K^y and L^y (which, of course, implicitly contain the temperature). 

Fig. 6 shows a small section of the Labyrinth tiling with the underlying grid (light) and 
the rapidity lines (dashed) with the associated parameters. Here, the rapidity lines are 
constructed parallel to the underlying grid in such way that they divide the edges of the 
underlying grid into two equal halves. Note that this gives a one-to-one correspondence 
between the intersections of two rapidity lines and the bonds of the tiling on which the 
intersection occurs. The checkerboard structure is apparent as the vertices of the Labyrinth 
lie in rectangles which constitute one of the two sublattices of the rectangular grid formed 
by the rapidity lines. 

It is not possible to give a full account of the actual solution in what follows. Instead, 
with special focus on the Labyrinth, we summarize the key points of the approach of Ref. 
P^] which the interested reader should have at hand for cross-reference. Though we can 



simplify several of the conditions (and actually give the results without derivation) we still 
depend on the explicit parametrization by elliptic functions which we will now outline. 

4.1 Elliptic parametrization and duality 

One can explicitly obtain parametrizations of the couplings in the integrable case for three 
different regimes ||36[, namely 



sinh(2ir,^) = F{ri-u', + Uy), sinh(2L,^) = F{u',-Uy) (4.1) 

where F{u) is a quotient of Jacobi's elliptic functions sn(-u), cn('u), and dn(-u) and where tj 
is determined by their half-period magnitudes K and K' , see Table 1. The rapidities u^ and 
u'^ {x G {a, &}) are a priori arbitrary complex numbers, but in order to obtain real (positive) 
coupling constants they have to be chosen appropriately, see Ref. |3^ for details. In Table 1, 
we also define several quantities which will be used in what follows. Here, k and k' denote 
the elliptic modulus and conjugate modulus, respectively, and Q{u) and 6i(u) are the usual 
Jacobian theta functions W^ ■ 
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Table 1: Parametrization in terms of elliptic functions 





regime I 


criticality 


regime II 


regime III 


n 


1 
k' 


1 


k' 


ik' 
k 


V 


K 


n 
2 


K 


K-iK' 


F{u) 


sn(w) 
cn(u) 


tan(M) 


cn{u) 


cii(uj 


r{u) 


dn(u) 
sn(ii)cn(M) 


2 


dn(u) 
sn(M)cn(M) 


cn(u) 
sn('u)dn(w) 


sin(2u) 


q{u) 


K'Q'M 

TT Ql{u) 





K' e'{u) 

TT e{u) 


K' Q'{u) 

TT Q{u) 


c 


n 





k^K' 

TT 


K' 

TT 



In all regimes, we require the elliptic moduli k, k' {k' = y/1 — k^) to satisfy the inequality 
< /c, /c' < 1. Hence, the regimes are characterized by the values of Vt (which plays the role 
of a "temperature- like" variable) given in Table 1 as follows: on has D? > 1^ < i7^ < 1 
and D? <Q for regimes I, II and III, respectively. The critical surface {D? = 1) separates 
the "low-temperature" regime I from the disordered "high-temperature" regimes II and III. 
At all other values of Q, the partition function turns out to be analytic. We have listed 
regime III for the sake of completeness because it might become interesting in the context 
of Lee- Yang zeros in the complex temperature plane. 

In the present context, we are mainly interested in regimes I and II as these lead to 
positive real couplings for the Ising model. In particular, this is the case if one chooses real 
rapidity differences u'^ — Uy within the range < u'^ — Uy < rj. Obviously, these two regimes 
are also related by duality (which inverts fi^) since the dual couplings (defined by Eqs. ( |3.22|) 
and ( p.23|) ) are obtained from 



sinh(2is:^ 



xyJ 



n~^F{r^-u', + u 



smh{2Lly) = n~' F{u', - Uy) 



(4.2) 



Note that in the parametrization given in Table 1 the critical case corresponds to the limits 
k ^ 0, k' —>■ 1, K —>■ 7r/2, K' —>■ oo, and q ^ where q = exp{—7cK'/K) is the elliptic nome. 
The above parametrization defines a subspace of couplings for which the Ising model is 
"Z-invariant" . The main ingredient in the derivation is, of course, the star-triangle relation 
P, ^. Although we do not want to go into details of the calculation, we still feel that a 
few words should be spent on the actual conditions on the couplings which are solved by the 
parametrization (^4.1|). 
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The restrictions imposed by the condition of integrabihty can be conveniently summarized 
as follows. To this end, it is convenient to introduce another involution ? , 



^xy J \ L^y 



(4.3) 



where Kr,,y and L^y are implicitly given by 

s\n\i{2K^y) smh{2K^y) = 1 , sinh(2ir^y) sinh(2L^y) = 1 . (4.4) 

This involution, together with the duality transformation * , generates the Abelian group 
Z2©Z2 (Klein's 4-group) which acts on the 8D coupling space IR^. (Note that ? and * 
coincided in Sec. 3.4 which was the reason for the "poor" self-duality condition). In the 
present context, integrabihty means that the equation 

^2 _ yr sinh(2M^) cosh(Mi + Ma + M3 + M4 - 2Mj) 
/=! sinh(2Mj) cosh(Mi + M2 + M3 + M4 - 2Mj) 

has to be satisfied for all vertices, where Mj (j = 1, . . . , 4) denote the couplings on the four 
surrounding edges of the vertex. In the Labyrinth, we have nine different vertices (counting 
orientations), cf. Fig. 4, which yield five different equations (due to symmetries). Note that 
the value of fl is the same for all vertices, hence one equation can be regarded as the definition 
of Q and the remaining four result in relations between the couplings. Using the notation of 
Eq. ( |3.24| ), four of the five conditions can actually be simplified to 

Saa = Sab = Sha = S^b = ^ (4.6) 

which is also apparent in the parametrization given in Table 1 as the function F satisfies 

F{u)F{r]-u) = n (4.7) 

in all three regimes. 

The one additional equation stems from the vertex which involves one long, one short 
and two medium length bonds, see Fig. 4. The equivalence of the equations obtained from 
the four different orientations follows from Eq. ( |4.(j| ). 

Note also that Eq. (|4.5| ) is automatically satisfied for vertex configurations obtained by 
including words with sequences . . .aa . . . which do not occur in the silver mean chain (|2.1| ). 
This means that the same conditions arise if one considers Ising models on tilings which 
are defined by different two-letter substitution rules in an analogous way. If one wants to 
consider models which involve only the four conditions of Eq. (|4.6| ), one is forced to use 
periodic sequences of just one letter in one direction as this is the only way to avoid the 
vertex which yields the additional equation. 

The conditions ( ^3| ) drastically restrict the possible couplings for the integrable cases. 
In the "isotropic" case Kaa = Laa = Ks, Kab = Lab = Kba = Lba = Km, Kbb = Lbb = Ki con- 
sidered in Sec. 3, the solvable subspace reduces to the isotropic Ising model on the square 
lattice, since from Eq. ( [4.6|) one obtains smh.{2Ks) = smh.{2Km) = sinh(2i^£) = Q^/^. 
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Of course, Eqs. ( [4 .51 ) and ( [4 .61) are entirely consistent with the results of Sec. 3 as the 
intersection of the solvable subspace with the self-dual surface ( |3.26| ) is just the critical 
solvable case with fi = 1. One should also realize that it is the additional equation with 
respect to Eq. ( ^^ that actually restricts the solvable critical manifold to a subspace of the 



self-dual manifold in our case. If one considers models where this equation is absent (which 
means that one uses a periodic sequence aaa... or bbb... in one direction of the grid on which 
the tiling is built), the two spaces coincide. It would be interesting to see explicitly how the 
extra condition in the generic case affects the critical behaviour of the system. 

4.2 Partition function and magnetization 

Provided that Q and all couplings are real and positive, we can follow the techniques outlined 
in |[l^, ^, |18|, |19| to calculate the partition function in the thermodynamic limit. Here, the 
result (normalized per site) reads 

Z := jimZy^" = (4fi)i/^expQ/i,U + ^/irn(a + ^6a) + ^M66) (4.8) 

where 

^xy = 0« - Uy) + (j){ri - < + Uy) (4.9) 

and (cf. Table 1) 

1 r / X 



(j){u) = -cu + / 7777 + 5' (a^) r{x)dx . (4.10) 

2 JO \2K 



The local spontaneous magnetization {a) is given by |3 

^ ^ \ if f]2 < 1 ^ ^ 

for any single Ising spin of the lattice - wherefore we suppressed the site index in ( [4.11| ). 
One can also easily find expressions for the correlations of two Ising spins that belong to 
the same plaquette of the Labyrinth tiling |3^. These, of course, depend on the particular 



couplings along the edges of that plaquette. 

The above feature is certainly somewhat surprising as one would expect the one-point 
function to depend on the local neighbourhood. That this is not the case for the magneti- 
zation is clearly due to the conditions (|4.6| ) on the couplings. It is to be expected that one 
would, for instance, observe a smaller magnetization at the vertices which are surrounded 
by long bonds only if one chooses Kbb and L^b to be small simultaneously (compared to the 
other couplings) which necessarily violates Eq. ( [4.61) . In this sense, the integrable subspace 
is possibly not representative and we presently cannot exclude the possibility of more than 
one phase transition in the general case where the corresponding order parameters would be 
the magnetization of certain subsets of the Labyrinth vertices. 
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5 Concluding remarks 

A zero-field Ising model on the Labyrinth tiling has been investigated. We used duality 
arguments to obtain information about the location of critical points and determined the 
self-dual couplings. Also, the subclass of exactly solvable models was considered in detail, 



following closely the discussion of the checkerboard Ising model in Ref. |^. In this case, the 
self-dual couplings also define the critical surface, the transition being of the Onsager type 
as in the periodic model. 

The special model investigated here is, of course, just an example. It can be generalized 
to tilings obtained by an analogous procedure from arbitrary two-letter substitution rules, 
provided one uses the correct periodic approximants. In this way, one can study Ising systems 
where the fluctuations of the underlying tiling behave in different ways. This is of interest 
because it has recently been shown in several closely related examples that these fluctuations 



actually determine the critical behaviour of the Ising model, see e.g. [^, [T^, [Tj, |T^. 

Of course, the investigation of the solvable cases can only be the first step in the investi- 
gation of the critical properties of Ising models on aperiodic tilings, and we hope to report 
on further results soon. The exactly solvable Z-invariant case is arguably too restrictive to 
be assumed representative as it, for instance, forces the magnetization to be independent of 
the location, which will certainly not be the case for some arbitrarily chosen couplings. 

Finally, one can extend the analysis to other systems like the more general eight-vertex 
model in the spirit of [|r^, R, T^. Here, it is now relatively clear what happens in the solvable 
cases which are, by no means, restricted to de Bruijn grids, cf. [^]. However, any progress 
beyond local solvability would help understanding the hierarchy of critical phenomena be- 
tween periodic and random. 
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Figure 1: The Labyrinth with underlying grid and boundary conditions 
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Figure 2: The Labyrinth and its dual 
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Figure 3: The tiling and its inflation 
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Figure 4: The tiles and their dual vertex configurations 
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Figure 5: Assignment of couplings K^y and L^y {xy G {aa, ab, ba, bb}) 
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Figure 6: The tiling with rapidity hnes and underlying grid 
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